Quantum Trajectories for Realistic Photodetection I: General Formalism 
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Quantum trajectories describe the stochastic evolution of an open quantum system conditioned on 
continuous monitoring of its output, such as by an ideal photodetector. In practice an experimenter 
has access to an output filtered through various electronic devices, rather than the microscopic 
states of the detector. This introduces several imperfections into the measurement process, of 
which only inefficiency has previously been incorporated into quantum trajectory theory. However, 
all electronic devices have finite band widths, and the consequent delay in conveying the output 
signal to the observer implies that the evolution of the conditional state of the quantum system 
must be non-Markovian. We present a general method of describing this evolution and apply it to 
avalanche photodiodes (APDs) and to photoreceivers. We include the effects of efficiency, dead time, 
bandwidth, electronic noise, and dark counts. The essential idea is to treat the quantum system and 
classical detector jointly, and to average over the latter to obtain the conditional quantum state. 
The significance of our theory is that quantum trajectories for realistic detection are necessary for 
sophisticated approaches to quantum feedback, and our approach could be applied in many areas 
of physics. 

PACS numbers: 03.65.Yz, 03.65.Ta, 42.50.Lc, 42.50.Ar 



I. INTRODUCTION 
A. Measurement in Quantum Mechanics 

To obtain information about a system of interest a 
measurement has to be made. In experiments that probe 
the quantum nature of our world, the system itself is, 
in general, necessarily affected by the act of measure- 
ment. For experiments involving continuous monitoring, 
traditional (projective) measurement theory as axioma- 
tized by von Neumann jjj is inadequate as a description. 
Such measurements would prevent any interesting evolu- 
tion occurring, because of the quantum Zeno effect [Q. 
To avoid this, a generalized quantum theory of measure- 
ment must be used |3|, [| . Generalized quantum measure- 
ments can be derived by considering projective measure- 
ments on a second system interacting with the system 
of interest. Often the second system has a clear phys- 
ical interpretation. For example, in quantum optics a 
bath (continuum) of electromagnetic field modes plays 
this role and the quantum system is said to be open. 

If the system is weakly coupled to its bath and the dy- 
namics are such that information concerning the system 
is spread throughout the many degrees of freedom of the 
bath, then a Markovian evolution equation for the sys- 
tem alone may be obtained. This equation is known as a 
master equation (ME). If, in addition, the bath is being 
measured then a stochastic master equation (SME) for 
the conditioned state of the system can be found. This 
is termed a quantum trajectory ||. The relation of the 
ME to quantum trajectories is that the former results 
from summing all possible trajectories, weighted by their 
probability of occurrence. Thus, the ME represents av- 
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erage evolution and gives the evolution of the state of 
the system when the environment is not monitored. A 
quantum trajectory (which is said to 'unravel' the master 
equation) gives the evolution of the state of the system 
conditioned on the results of monitoring the environment. 
A quantum trajectory is necessarily stochastic, and can 
be jumpy |, |, g or diffusive §, | |, @§ @, de- 
pending on the nature of the measurement. 



B. Realistic Quantum Trajectories 

In experiments performed in the laboratory, a detec- 
tor mediates between the detector input (for example, 
the electromagnetic field emitted by a fluorescent quan- 
tum system) and the experimenter. Realistic detectors 
are not perfect. Information is lost in the conversion 
of the quantum field to a signal that the observer can 
use. This loss may occur at the front end of the detector 
(characterized by an inefficiency), in the circuit contain- 
ing the detector (described by a response time and elec- 
tronic noise) or at the circuit output (electronic output 
noise). In order to obtain a true quantum trajectory for 
the experiment, the observer must condition the state 
of the quantum system on results that are available in 
the laboratory rather than on the microscopic absorption 
events considered previously in quantum trajectories. 

Realistic quantum trajectories of this sort cannot be 
encompassed within the standard formalism of Refs. 

except for inefficiency |l3|]. This is because the elec- 
trical circuit, which filters the measured output of the 
quantum system, causes a statistical delay in the observer 
obtaining information about the system. Thus the con- 
ditioned evolution of the quantum system alone is non- 
Markovian. A method for treating this was first proposed 
by us and Mabuchi |I3], in which the quantum system 
is embedded within a supersystem that also contains the 
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state of the detector. If the set of (classical) detector 
states is S, then the supersystem is described by the set 
{p s : s e §}. Here Tr[,o s ] is the probability that the ap- 
paratus is in state s, and p s /Tr[p s ] is the system state 
given this event. The supersystem state, {p s }, obeys a 
Markovian evolution equation. The detector states are 
classical in the sense that it is not necessary to consider 
superpositions between them. 

In this paper we consider in detail two sorts of real- 
istic detectors, an avalanche photodiode in geiger mode 
(a photon counter) and a photoreceiver (for homodyne 
detection). The detector states in the former case are 
discrete (we use three states) and in the latter contin- 
uous (although discretization is necessary for numerical 
simulation). For the photon counter, this model allows us 
to include the effects of detection efficiency, dark counts, 
response (or rise) time and dead time. For homodyne 
detection using a photoreceiver we investigate detector 
efficiency, output Johnson noise (which adds onto the 
measurement result linearly) and detector bandwidth. 
(Of course, this bandwidth relates not to the range of 
frequencies (colours) of light to which the detector will 
respond, but to the characteristic time it takes for the ab- 
sorption of light by the detector to generate a response 
in the detector output.) 

Before passing onto our method of derivation of real- 
istic quantum trajectories, it is worth mentioning that a 
three-level model for a photon counter including many of 
the above effects has been considered by Gardiner p5[ . 
However this theory is based upon quantum Langevin 
equations rather than quantum trajectories, and does not 
allow one to determine the system state conditioned on 
the experimenter's knowledge. 
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FIG. 1: Schematic diagram for the method that we use to 
derive realistic quantum trajectories. It should be noted that 
some of the steps shown above can be performed in a dif- 
ferent order. This flexibility is useful for keeping the deriva- 
tion of realistic quantum trajectories as simple as possible. 
The following abbreviations have been used: stochastic mas- 
ter equation (SME), stochastic differential equation (SDE), 
stochastic differential Chapman-Kolmogorov equation (SD- 
CKE) , Kushner-Stratonovich equation (KSE) and superoper- 
ator Kushner-Stratonovich equation (SKSE). Hexagons, rect- 
angles and circles indicate bodies of knowledge, equations and 
procedures, respectively. See text for further details. 



C. Derivation of Realistic Quantum Trajectories 

The derivation of realistic quantum trajectories con- 
sists of a number of well defined steps, that are schemat- 
ically illustrated in Fig. [|. These can be understood by 
considering the subsystems that make up the supersys- 
tem, namely the quantum system and the detector. 

The basic philosophy that we adopt at the start of the 
derivation is that in principle all of the microscopic states 
of the detector could be monitored, without altering the 
quantum trajectory of the system. Essentially this is 
because the detector states consist of distinct configura- 
tions of an enormous number of particles. These particles 
are constantly interacting with their environment, thus 
making coherences between detector states unimportant. 
Thus in principle it is possible for an experimenter to 
have sufficient information to unravel the evolution of 
the system by a conventional stochastic master equation 
(SME). In practice, an experimenter has access to far 
less knowledge. This loss of knowledge is due to classical 
uncertainties introduced by the detector. 

To describe the detector we consider some property of 
it that can be related to the output that the observer 



monitors. For example, in the case of homodyne detec- 
tion the charge across one of the capacitors in the cir- 
cuit represents the detector state. By considering all of 
the dynamic influences upon this circuit variable we for- 
mulate an evolution equation for it that is stochastic in 
general. Since finally we will not know the stochastic 
inputs, we transform this into a stochastic differential 
Chapman-Kolmogorov equation (SDCKE) that follows 
the probability of the various states being occupied. We 
use this terminology to indicate a generalization of the 
deterministic differential Chapman-Kolmogorov equation 
discussed in O]. 

The realistic measurement results available to the ob- 
server give information about the detector state, without 
necessarily revealing it fully. Based on the result and the 
a priori (initial) probability distribution for the detector 
state, the best estimate of the updated state (the a pos- 
teriori distribution) is found. This is done via a Bayesian 
analysis Q . The resulting equation for the the probabil- 
ity distribution is, for the case of Gaussian white noise, 
usually referred to as a Kushner-Stratonovich equation 
(KSE) |jl8). For convenience we will use this terminology 
even in the case of point-process noise. 
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Note that normally a KSE would also contain that un- 
conditioned (deterministic) evolution due to the system 
dynamics. Since these dynamics are stochastic at this 
stage of the derivation (they are contained in the SD- 
CKE), we do not include them in the KSE for the de- 
tector state. Note also that the SME, despite its name 
(introduced in Ref. ]l3|]), is conceptually closer to the 
KSE than to the SDCKE: the stochasticity in both the 
KSE and the SME is solely due to measurement noise, 
and the stochastic terms are nonlinear in both cases. In 
fact, the SME can justifiably be called a quantum KSE 
& 

The evolution of the SME, the SDCKE and the KSE 
are combined to give the joint stochastic equations gov- 
erning the evolution of the supersystem. Up until this 
stage all the stochastic influences on the supersystem are 
still present. To link the evolution of the system with 
that of the detector this is necessary. Now, however, we 
average over processes that a realistic observer cannot 
monitor. This leaves only the stochasticity associated 
with the measurement and our derivation of the realistic 
quantum trajectory is then complete. We term this final 
stochastic trajectory for the supersystem a superoper- 
ator Kushner-Stratonovich equation (SKSE). Obviously 
the reader will gain much greater insight into this process 
by following the examples provided in the later sections. 



D. Importance of Realistic Quantum Trajectories 

Thus far, the main practical utility of quantum tra- 
jectory theory has been in improving the computational 
efficiency of simulations used to compare models with 
experimental data. For this purpose, the introduction 
of realistic quantum trajectories could only hinder sim- 
ulations. However, quantum trajectory theory is now 
gaining increasing importance as the quantum gener- 
alization of Kalman filtering, which provides essential 
signal-processing methods in classical estimation, com- 
munication, and control engineering. Quantum trajec- 
tory theory should in principle play the same pivotal 
role for emerging quantum analogs of these technologies 
p0| , pl| p2[ . Before this can happen it is essential that 
the theory be extended to account for the imperfections 
of realistic measurement devices. Non-ideal detector dy- 
namics can dramatically affect the proper inference from 
measured signals to the conditional quantum state of an 
observed system, as we will show. 



state. A discussion of non-realistic (that is, ideal) quan- 
tum trajectories follows this, which allows the SMEs for 
the quantum system to be found. 

In Sees. Ill and |iy|, realistic quantum trajectories are 
derived for photon counters and photoreceivers respec- 
tively. Both of these derivations are based on the pro- 
cedure outlined in Sec. IC. For the photon counter, we 
finally idealize the measurement process to show that our 
equations reduce to those associated with perfect detec- 
tion (traditional quantum trajectories). An important 
and non-trivial consideration in the case of the photore- 
ceiver is the bandwidth of the device. As the photore- 
ceiver circuit adds electronic noise, a calculation of this 
bandwidth based on just the resistance-capacitance (i?C) 
time constant of the circuit does not yield the correct 
answer. As a final point we consider the effect of letting 
the RC time constant go to zero but still maintaining the 
presence of the electronic noise which obscures the input 
signal. We conclude with a summary of our results and 
discussion of future directions in Sec. [v|. A numerical 
exploration of our theory is given in the following paper 



II. STOCHASTIC EQUATIONS 

As discussed in Sec. IC, there are three types of 
stochastic equations that go into making a realistic quan- 
tum trajectory. In this section we examine the general 
theory for each. 



A. Stochastic Differential Chapman-Kolmogorov 
Equation 

A typical Langevin equation or stochastic differential 
equation (SDE) for a classical random variable X can be 
written as 



dX = a(X)dt + VDdW{t) +edN{t). (2.1) 

Here a(X) is an arbitrary function of X, and dW(t) 
and dN(t) are stochastic increments. The point pro- 
cess dN(t) is either zero or one, and has an infinitesimal 
mean which is a positive function of X , say g(X)dt. The 
Wiener increment dW(t) is related to Gaussian white 
noise £(i) by dW(t) = ^{t)dt. Gaussian white noise has 
statistics of |16i 



E. Overview of this Paper 

As stated above, the main thrust of this paper is to de- 
rive in detail the realistic quantum trajectories for pho- 
ton counting and homodyne photodetection introduced 
in Ref. H|. In Sec. [n], we first briefly give the theory of 
classical trajectories and measurement, which is impor- 
tant for obtaining SDCKEs and KSEs for the detector 



EK(t)] = o, = s(t f). 



(2.2) 



Here the roman font E denotes exp ecta tion value. By 
writing the SDE in the form of Eq. (2.1) we imply that 
the equation is to be interpreted in the Ito Q or explicit 
Ipif sense. This means, for example, that the expected 
value of dX may be found by replacing dN and dW by 
their expected value. This is not true for the implicit 



or Stratonovich |16 form. By the same token, the chain 
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rule of standard calculus does not apply in general to the 
Ito form, whereas it does for the Stratonovich form. 

Langevin equations are appropriate for situations in 
which the stochastic increment is known, thus allowing 
X to be tracked. Often it is relatively easy to turn the 
physics of a problem into a Langevin equation. An ex- 
ample of this is a noisy electrical circuit. The analysis of 
such circuits is well understood and is based on simple 
principles such as Kirchhoff's laws. 

In this paper, Langevin equations are first found on the 
presumption that the stochastic increment is known. In 
reality, however, this is not the case and we cannot track 
the variable X perfectly. For this reason it is necessary to 
turn the Langevin equation into an equation describing 
the evolution of the probability distribution for X, which 
is written as P(x). Note that the transformations that we 
perform actually maintain the stochastic nature of this 
evolution equation, as the averaging over unobserved pro- 
cesses is done at a later stage. This is necessary because 
the noise in the equation for P(x) may also appear in 
other parts of the overall derivation. Before this final av- 
eraging is done, the stochastic equation for P(x), would 
preserve a Dirac <5-function as a solution. We now do this 



for Eq. (2.1), as an aid the reader. 



If X is known at time t, then at that time P(x) = 
S(x — X). This then evolves according to 

dS(x -X) = 5{x-X- dX) - S(x - X) (2.3) 
= 5[x-X- a(X)dt - VZ> dW(t) - edN(t)} 
-5(x-X). (2.4) 

We can treat the diffusion and jump noise in a u nified 
fashion by constructing a Taylor series of Eq. (^4) to all 
orders of the first term on the right-hand side (RHS): 



dS(x - X) 



exp 



d_ 

dx 



edN(t) 



-a{X)dt-VD dW(t) 



l)S(x-X). (2.5) 



We now expand the exponential to all orders and use the 
Ito rules dN 2 = dN, dW 2 — dt, with all other products 
being zero. This gives 



dS(x 



' dx 



a(X)dt + VDdW(t) 



D 92 M \ u 
+ Jdx^ d T {x 

+ dN(t) [S(x-X- 



X) 



5(x-X)] 



(2.6) 



Using the delta function to change a(X) to a(x) and then 
averaging over the random variable X gives the Ito equa- 
tion for the probability P(x) = E[5(x — X)], 



d 

dx 
x P(x 



i D d 2 
-a(x)dt - VDdW(t) + —-^dt 
J 2 dx z 

)+dN(t)[P(x-e)-P(x)]. (2.7) 



The first (dt) term is deterministic drift, the second (dW) 
is stochastic drift, the third (dt) is diffusion and the last 
(dN) is the jump term. We refer to this equation as 
a stochastic differential Chapman-Kolmogorov equation 
(SDCKE). If there were no jump term, then we would 
have a stochastic Fokker-Planck equation. 

If averages over dW(t) and dN(t) are taken then we 
have a standard (deterministic) differential Chapman- 
Kolmogorov equation 



d . . D d 2 1 . , 
■g(x-e)P(x-e)-g(x)P(x), (2. 



where g(x) enters from the definition ~E[dN(t)f(X)] 
E[g(X)f(X)}. 



B. Kushner- Stratonovich Equation 

From the point of view of measurement, it useful to 
consider the state of a classical system to be defined by 
the probability distribution on state space. This is the 
space of all relevant physical quantities, called state vari- 
ables. In this paper the classical system will be a detector 
and the associated electrical circuit. Only one variable is 
necessary to describe it so we can think of the distribu- 
tion P(x) that was introduced in the previous section as 
representing the state. 

As the state of the system summarizes the observer's 
knowledge it is logical that the state will change on the 
basis of measurement upon the system. From the a priori 
state P(x) of the system and the result r of measurement, 
the a posteriori state P(x\r) can be found using Bayesian 
inference |l7j 



P(x\r) 



P(r\x)P(x) 
P{r) 



(2.9) 



Here P(x\r) reads as 'the probability of x given r\ The 
division by P(r) is for normalization purposes, where 



P(r) = / dxP(r\x)P(x). 



(2.10) 



Here P(r\x) represents the probability of obtaining re- 
sult r, g iven that the state is x. The expression in 
Eq. ( |2.9[ ) will be repeatedly used in this paper to obtain 
the Kushner-Stratonovich equations (KSEs) for the clas- 
sical system. This is simply the continuous-in-time limit 
of Bayes theorem. Rather than provide the general form 
of the KSE here, the reader will be guided through the 
derivations of the KSEs (as they arise) from fundamental 
considerations. 

The unconditioned system state which we would use 
if the measurement had been performed and ignored is 
given by 



P'(x) = / drP(x\r)P(r), 



(2.11) 
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where we have assumed that the results take on contin- 
uous values. The prime is included to remind the reader 
that this is an a posteriori state. In the above we have 
assumed that there is no measurement back-action on 
the system - that is, the measurement only affects the 
observer's knowledg e and not the dynamics of the sys- 
tem itself - so Eq. (2.11) reduces to P(x) as one would 
expect. Of course, it is possible to consider classical mea- 
surements for which there is a measurement back-action. 
This is also the case for all quantum measurements be- 
cause of the Heisenberg uncertainty principle. 



C. Stochastic Master Equation 

As noted in Sec. I C, the SME is a quantum version of 
the KSE, including the evolution unconnected with the 
measurement as well. The SME can be derived within the 
standard formalism of generalized quantum measurement 
(for simple introductions, see for example Refs. |^5| and 
p6[ ). This formalism is simply the quantum analogue 
of Bayes' theorem, generalized for measurement back- 
action, as shown recently in Ref. |27|. In this section 
we will not review the derivation of the SME from the 
standard formalism, but will rather just quote the final 
results. 



where the superoperator T> is defined, for arbitrary B, as 



1. Jumpy Trajectories 

In many examples of continuous measurement on the 
output system, the result can only take on two possible 
values. This is the case in ideal direct photodetection 
in which the electromagnetic field is monitored with an 
ideal detection device that responds to the presence of 
photons. In any infinitesimal time interval a photon ei- 
ther is or is not detected. Such measurements give rise 
to jumps in the conditional system state because the rare 
event of a photon detection conveys to the observer a fi- 
nite amount of information about quantum system. For 
example, if the fluorescent system were a TLA then a 
photon detection implies that the TLA now resides in 
the ground state. A null result still causes a non-unitary 
change to the system state, but it is of infinitesimal mag- 
nitude. That is not to say that it is unimportant, how- 
ever, as over a finite length of time between jumps, it will 
also cause a finite change in the system state. 

Let the fluorescent system have lowering operator c, 
and assume that the fluorescence can be formed into a 
beam of light. Also let an optical local oscillator (that is, 
an effectively classical light field from a laser) be added 
to the fluorescence before it is detected. This can be done 
using a highly transmitting beam-splitter. Let the local 
oscillator (LO) amplitude added to the fluorescent beam 
be p. If we measure time in units such that the photon 
flux from the system is (c^c), then the unobserved system 
obeys the (deterministic) master equation 

p(t) = -i[H, p(t)} + V[c]p(t) = Cp, (2.12) 



V[B]p = J[B]p-A[B]p 

= BpB^ - i (B^Bp 



pB^B) 



(2.13) 
(2.14) 



where the two terms in the last line define J and A 
respectively. This master equation is of the Lindblad 
form |28| , where we have assumed for simplicity that the 
evolution apart from the radiation damping term (T>[c]p) 
can be described by a Hamiltonian H. 

To give an explicit description of the conditional sys- 
tem state evolution under ideal photon detection, it is 
useful to define the point process, dN(t), that is equal 
to 1 if there is a detection in the time interval [t, t + dt) 
and otherwise. Allowing for a detection efficiency rj, 
the statistics of dN(t) are defined by 



dN(t) = dN{t) 2 



E[dN(t)] = 77Tr[(ct + /j*)(c + fi)p N ]dt 



(2.15) 
(2.16) 



The subscript N on the system state is to indicate that 
it is conditioned on the detection events. The evolution 
of the system state in terms of dN(t) is given by |2~|| 

dpN = — dtH[iH + ^rjc'c + rjp*c + \rj\p\ 2 ]pN 
+ dN(t)g[^Tj{c + p)} PN 
+ dt(l-rj)V[c]p N . (2.17) 

This equation is a stochastic master equation (SME). 
Note that we have dropped the time argument of the 
conditioned system state, but maintained that of the in- 
crement, dN(t). This is done for convenience and, also, to 
emphasize the stochastic nature of the detection events. 
Time arguments are only included in this paper when 
they will aid the reader's understanding. The non-linear 
superoperators, Q and Tl are defined by 



G[A]B = 



J[A]B 



- B. 



(2.18) 



Tc[J[A]B] 

H[A]B = AB + BA^ -Tr[AB + BA^]B. (2.19) 



It is possible to return to t he ME (Eq. ( |2.12|) ) by re- 
placi ng dN (t) in Eq. ( [2.17| ) by its expectation value 
(Eq. (gig)). 



2. Diffusive Trajectories 

In the limit of p — > oo the rate of quantum jumps goes 
to oo but the effect on the system of each one goes to 
zero. This leads to quantum trajectories of a diffusive 
nature, where the system evolves continuously but non- 
differentiably in time. In this limit the photocounts give 
way to a photocurrent, J. Including a detection efficiency 
77, and for a particular normalization, it is given by 



Jdt = ri (ce _i * + cV*) dt + y/fjdW(t), 



(2.20) 



where dW(t) = £(t)dt represents Gaussian white noise 
p6| having the properties given in Eqs. (|2.2[)-(|2~2]). Here 
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REALISTIC PHOTON COUNTER 





FIG. 2: Schematic diagram of the model used in this paper 
for a realistic photon counter monitoring a TLA placed in an 
optical cavity, the output of which is combined with that of a 
weak local oscillator (LO) at a low reflectivity beam splitter 
(BS) before direct detection. This will be the system on which 
investigation will focus in future numerical work. In this di- 
agram we have used a BS of efficiency r\ followed by an ideal 
photon counter to represent the quantum efficiency of the re- 
alistic detector. The three states of the detector 0,1,2 are 
explained in the text. Single arrow heads are used for Pois- 
son processes. The photocurrent produced by the detector I 
will consist of spikes at the time of avalanches {dN = 1). 



( pop 



corre- 
while 



<J> = arg(/i), and putting $ = in Eq 
sponds to measurement of the x quadrature c + c 
$ = 7r/2 corresponds to measurement of the y quadra- 
ture —i(c— c'). The system state evolution conditioned 
on this photocurrent is || |l3| 

dp j = dtCpj + {J(t)dt - E[J(t)]dt} H[e~ iiSf c]pj, (2.21) 

where again the subscript on pj indicates that it is con- 
ditioned on the recording of the current J. 

III. REALISTIC PHOTON COUNTING 



FIG. 3: Diagrams (A) and (B) show electron energy band dia- 
grams for a p-n junction operated under a beyond breakdown 
reverse bias. In diagram (A) an incident photon causes the 
creation of an electron (filled circles) - hole (circles) pair that 
begins to separate due to the applied electric field. The APD 
would be in state 1. In (B) the electron and hole collide with 
other charges leading to an avalanche. Whether the avalanche 
has reached threshold or not determines if diagram (B) cor- 
responds to the detector being in state 2 or 1, respectively. 



field, leads to an avalanche forming. At this stage (be- 
fore the avalanche has been detected) the APD is still 
in state 1. The avalanche continues to build until the 
current reaches some threshold value and a detection is 
registered, thus changing the state of the APD to 2. Due 
to the stochastic manner in which the avalanche spreads, 
the dwell time in state 1 is not deterministic. We model 
the 1 to 2 transition as a Poissonian process with rate 7 r . 
We call the mean dwell time, 7" , the 'response time'. 
The avalanche, once detected, is arrested by the appli- 
cation of a negative-going voltage pulse that tempor arily 
brings the bias voltage below the breakdown value |pl| . 
This results in a hxed 'dead time', rjd, during which the 
APD cannot detect photons, after which it is restored to 
state 0. 



An avalanche photodiode (APD) operating in Geiger 
mode produces a macroscopic current pulse in response 
to an incident photon. This allows the observer to detect 
the presence of photons and is the reason the device is 
also known as a photon counter. The size of the current 
pulse is large compared to the sources of noise within the 
detection circuitry, which allows the threshold current 
value for a detection to be set well above this noise. 

The APD essentially consists of a p-n junction oper- 
ated under a reverse bias greater than the breakdown 
voltage |H, [50|, . Under these conditions the diode can 
be described by just three classical states |30| (see Fig. ^ 
and Fig. ||). This will enable a mathematical model of 
the APD to be simply constructed. 

The first state (0) is a stable low-current state in which 
there are no charge carriers in the depletion region of the 
junction. The transition from to the "unstable" in- 
termediate state (1) takes place when an electron-hole 
(e~-h + ) pair is created in the depletion region by an 
incident photon (with quantum efficiency 77) or by ther- 
mally initiated 'dark counts' occurring at a rate 7dk- Fur- 
ther impact ionization, under the influence of the electric 



A. Physical Explanation of Parameters 

Now that our basic model for the APD has been in- 
troduced, a further physical explanation is given for the 
quantum efficiency, response time, the dead time and 
dark counts. 

The quantum efficiency, 77, is largely determined by 
two considerations. These are the fraction of incident 
photons that are absorbed in the active region of the 
junction and the fraction of e~-h + pairs that lead to an 
avalanche (30| . The first is determined by the absorption 
co-efficient of the material at the wavelength of interest 
and, also, the thickness of the depletion region. If the ab- 
sorption co-efficient is too large then the photons will be 
absorbed near the surface, where there exists only a weak 
biasing electric field. Too small a co-efficient allows the 
light to pass straight through. The second consideration 
is a strong function of the overvoltage (beyond break- 
down), with higher fields ensuring that the creation of a 
charged pair leads to an avalanche. Unfortunately, this 
also leads to approximately exponentially increased dark 
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counts p2| . 

The response rate, 7 n is also a strong function of 
the overvoltage due to its effect on the rate at which 
charge carriers move through the junction. Other impor- 
tant considerations are the dimensions of the junction, 
which alter the time taken for the avalanche to spread 
throughout. Two important mechanisms that facilitate 
the spreading are the direct electrical diffusion of carriers 
into adjacent regions of the junction and the emission of 
photons from hot carriers that are absorbed in a different 
region ||§. 

The dead time, Tdd , is generally quite long as all charge 
carriers must be swept out of the active region before the 
overvoltage can be re-applied without another avalanche 
resulting. The application of a negative-going voltage 
pulse (active quenching) is more effective than letting a 
resistance in series with the junction reduce the reverse 
bias below breakdown as the current increases (passive 
quenching). The active quenching of the avalanche can 
be initiated as soon as soon as the current reaches a 
threshold value. Minimizing the circuit noise allows the 
threshold to be lowered thus reducing the response time 
and, also, the duration of the voltage pulse. 

The dark count rate, 7dk, can be reduced to insignifi- 
cance with the use of a high quality Silicon APD f{0|| . 
However, such devices are not appropriate for detec- 
tion at all wavelengths. In other materials the exis- 
tence of free-carrier trapping centers and thermally gen- 
erated carriers creates a considerable problem. During an 
avalanche, carriers can become trapped in a region of the 
junction, only to be released after the above breakdown 
voltage has been restored (causing an afterpulse). This 
leads to a dark count which is signal-dependent as the 
number of trapping centers generally increases with the 
frequency of avalanches. Cooling will reduce the thermal 
dark count. These complexities are not introduced into 
the APD model in this paper. 



B. Realistic Quantum Trajectories 

Our aim is to derive the quantum trajectories for the 
quantum system (the source of the light entering the 
APD) conditioned on the observation of an avalanche. 
In order that the derived equations be applicable to dif- 
ferent measurement schemes, we consider the output field 
of the system of interest to be combined at a low reflec- 
tivity beam splitter with that of a weak local oscillator 
(LO), having transmitted amplitude fi. If the LO is not 
weak, and contributes a large photon flux, then a photon 
counter is not the appropriate type of detector. On the 
other hand, /i can be set to zero to obtain direct detection 
without problems. 

In Fig. [l] the general method of deriving realistic quan- 
tum trajectories was given. The point that some of the 
steps in the method are interchangeable was made in the 
caption to that figure. For the APD, we will in fact per- 
form the conditioning upon measurement after the un- 



observed processes are averaged over. This will keep the 
equations that we manipulate during the derivation as 
simple as possible. 



1. Stochastic Master Equation 

As described in the introduction, the detection device 
can be treated in a classical manner. This implies that 
the presence of an observer monitoring the detailed be- 
haviour of the charge-carriers in the p-n junction would 
not alter the detection results. The hypothetical ob- 
server could also, in principle, deduce the incidence of 
photons that were absorbed outside the active region of 
the junction. In this way, a list of the times at which 
photons arrived could be compiled. Of course, such sen- 
sitive measurements are currently out of reach due to 
technological reasons, but the point is that it is consistent 
with quantum measurement theory to assume a quantum 
jump unraveling of the ME. In fact, we even assume at 
this stage that the measurement can have unit efficiency, 
as it is conceptually more satisfying to explicitly include 
inefficiency in our description of the APD. Thus the fic- 
titio us co nditional dynamics of the system are those of 
Eq. fl2.17| ), with n = I. 



2. Stochastic Differential Chapman-Kolmogorov Equation 

We can now formulate equations of motion for the 
probability of each of the three classical states of the 
APD being occupied. Because the state space is discrete, 
averaging over the stochastic processes would necessarily 
yield a classical master equation jl6| and we could refer to 
the equations we derive in this section as a stochastic clas- 
sical master equation. However, to solidify for the reader 
the processes involved in forming a realistic quantum tra- 
jectory (described in Sec. IC), we refer to the stochastic 
probability equations as SDCKEs. The stochastic equa- 
tion would apply if all transitions were monitored, and 
that when an electron-hole pair is created a distinction 
can be drawn between those created via photon absorp- 
tion and those via thermal fluctuations (dark counts). 
Because of the classical nature of the APD dynamics, 
this is logically possible. Eventually we will average over 
processes that a realistic observer cannot monitor. 

There are a number of processes that contribute to the 
evolution of the occupation probabilities. Here, we will 
treat dark counts, avalanches, APD resetting and inci- 
dent photons from the monitored quantum system. Since 
all transitions are being monitored by a fictitious observer 
and we assume, for simplicity, that the initial state of the 
APD is known the evolution of the probabilities is de- 
scribed by jump terms only. If we use <i/Vdk(i) = 0, 1 to 
define the occurrence of transitions due to dark counts, 
dN'(t) to indicate transitions due to photons from the 
quantum system of interest and dAf'(t) for avalanches 
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then 

dP = -dN dk (t)-dN'(t)+djV'(t-T dd ) (3.1) 

dPi = dN dk (t)+dN'(t)-dAf'(t) (3.2) 

dP 2 = dN'{t)-dN'{t-T dd ). (3.3) 

From these SDCKEs it can be seen that photons and 
dark counts, which lead to the creation of an electron-hole 
pair, take the detector from state to state 1. Avalanches 
cause the transition 1 — > 2 and the resetting of the de- 
tector, with the use of the delayed time t — r ddl returns 
it to the ready state 0. 

The use of the prime on dN'if) indicates that this point 
process is associated with, but not the same as, dN(t) 
which appears in Eq. ( [2.17 ). In fact dN'(t) is equal to 
dN(t), except for having some of the ones removed. This 
leads to 



dN{t)dN'{t) = dN'(t). 
Their expectation values are related by 
E[dN'(t)] = r)P E[dN(t)]. 



(3.4) 



(3.5) 



This reflects the two ways in which it is possible to have 
a "missed detection" [dN(t) = 1, dN'(t) = 0], namely, 
inefficiency or the APD's not being in the ready state. 

The use of the prime on dj\f'(t) is to distinguish an 
avalanche from an observed avalanche, denoted by dM (t). 
We note that such a distinction could be necessary if the 
efficiency of avalanche detection was less than unity. Al- 
though we do not include an avalanche detection ineffi- 
ciency here we feel that the prime may aid the reader in 
understanding the need for the conditioning of the equa- 
tions upon measurement of the presence of avalanches, 
which will occur later. The statistics of dAf'(t) are de- 
fined by its expectation value 



E[dJ\f'(i)] = ^Pidt. 



This comes from our assumption that the avalanche's 
reaching its threshold value is a Poisson process. Finally, 
the expectation value of dN dk (t) is 



E[dN dk (t)} = ldk P dt, 



(3.7) 



reflecting another Poisson assumption, here for the dark 
counts. 



3. Joint Stochastic Equations 

The evolution of the quantum system is linked to that 
of the detector through dN'(t). To obtain realistic quan- 
tum trajectories for the supersystem, which includes the 
quantum system and the classical detector states, we 
form the quantity 

Pi = PiPN, (3.8) 
where pi is the unnormalized state of the system given 
that the detector is in state i. The normalized system 
state is 



PO + Pi + P2- 



(3.9) 



It is only possible to assume that the quantum and de- 
tector states factorize if there are no correlations between 
them. This will be the case until we average over unob- 
served processes. Thus the evolution of pi can be found 
from 



Pi + dpi = (Pi + dPi)(pN + dpN), 



(3.10) 



with 
Eqs. 



the use of the increments in Eq. (2.17) and 
(3.1 )— (3.3) . Conditioning upon the measurement 



(3.6) 



of d Af(t) will be done as a late r st ep. Substitution into 
Eq. J3.10D with the use of Eq. (O) gives 



Po(t- 


-dt) 


= po + Po [dN(t)g[c + 


p] - dtH[iH +\c)c + p*c 


+ h\p\ 2 \]pN 








- [dN dk (t) + dN'(t) 


-dN'{t-T dd ) + dN'{t)Q[ 


c + p}] Pn 


(3.11) 


Pl(t- 


-dt) 


= Pl +P l [dN{t)g[c + 


p] - dtH[iH + \c)c + p*c 


+ \\p\ 2 \]pn 








+ [dN dk (t) + dN'(t) 


-dJ\f'(t) + dN'(t)g[c + n] 


PN 


(3.12) 


p-l{t- 


-dt) 


= p 2 + P 2 [dN(t)g[c + 


p] — dtTi[iH + c + p*c 


+ \\P\ 2 ]]PN 








+ [dJ\f'(t)-dJ\f'{t- 


- T dd )] pn, 




(3.13) 



remembering that all time arguments not stated are t. and the APD. 
The linking of the quantum system and the APD is man- 
ifest in the dN' (t)Q[c+ p]pN terms. If we were now to av- 
erage over the emissions from the quantum system a cor- 
relation would still develop between the quantum state 
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4- Average Over Unobserved Processes 



This leads to the following three results for the states 



In order to obtain trajectories that are based on infor- 
mation that is realistically available in the laboratory, we 
must now average over unobserved processes. These are 
the photon emissions and dark counts. Strictly speak- 
ing, we should also average over d!\['(t) since dM(t) was 
stated as specifying observed avalanches. However, we 
have not yet performed the conditioning due to the de- 
tection of the avalanches so these transitions will be left 
for the moment. Doing the averages before the measure- 
ment does not change the final result, but does make the 
equations more simple. 



Th e re levant averages are Eq. ( 2.16 ), E q. ( 3.5) and 
Eq. (O). Substitution of these into Eqs. (|3.11|)-(|3"T5|) 



gives 

p (t + dt) = p Q + dt (C - 7 d k - r]J[c + M])Po 

+ dJV'(t-T dd )p (3.14) 

pi(t + dt) = p 1 +dt£p 1 +dt(j d - k + i]J[c + p])po 

-dJV'{t)p (3.15) 

p 2 (t + dt) = p 2 +dtCp 2 +dN'{t)p- dN' (t-T dd )p. 

(3.16) 



Note that p is given by Eq. ( |3.9| ) and C is specified in 
Eq. (|2T2|). 



5. Kushner-Stratonovich Equation 

The probabilities assigned to the occupation of the 
three APD states by the realistic observer will undergo 
jumps in the event of an avalanche. If there is no 
avalanche then this also provides information to the ob- 
server about the probabilities. Let us first consider the 
effect of an avalanche. The new states will be given by 



Tr[p, 



rP 



i\dj\f(t)=l 



(3.17) 



where Pi\dsf(t)=i is the probability of the APD being in 
state i given that an avalanche was observed. These prob- 
abilities can be simply calculated using Bayes conditional 
probability theorem JL7| 



Pi 



P 



i\dj\f(t)=l 



dM(t) = l\i"i 



P 



<W(t) 



(3.18) 



Here, PdU(t)=i\% is the probability of observing an 
avalanche given that the APD is in state i, Pi is the prior 
occupation of state i and P<iM(t)=i is the probability of 
observing an avalanche. The relevant probabilities are 

PdM(t)=i\o = PdM(t)=i\2 = (3.19) 
P<W(t)=i\i = Jrdt (3.20) 

%((]=! = IrPidt. (3.21) 



Po|djV(t)=i — 



Pi 



Pi\mt)=i ~ Tr[pi] 

P2\dM(t) = l = 0. 



(3.22) 
(3.23) 
(3.24) 



If there is no observed avalanche then we need to con- 
sider probabilities of the form Pi\dM(t)=Q- This can be 
calculated from 



P 



i\dN{t)=0 



PdM(t)=0\iPi 



P 



dM(t)=0 



(3.25) 



with the probabilities interpreted in the usual fashion. 
The relevant expressions are 

PdAf(t)=0\0 = PdAf(t)=0\2 = 1 (3.26) 

fW(t)=o|i - 1-Trdt (3.27) 
*W(t)=o = l-Tti\A- (3.28) 

The new states based on these probabilities are, to first 
order in dt, 

Pa\dN(t)=o = (1 + Ji-Pidt)p (3.29) 
Pi\dtf{t)=a = (1 - 7r<& + 7r-P]d*)/0i (3.30) 

P2\dN(t)=V = (l+JrPldt)p 2 . (3.31) 

We can summarize the avalanche and no avalanche con- 
ditioning with 

Po\dmt) = [l + lTPxdt-dN{t)}p (3.32) 
Pl\dM(t) = (1 - 7rd* + IrPidt) pi 

1 



-dAf(t) 



Tr[pi] 



1 Pi (3.33) 



p2\dN{t) = [1 + JiPidt - dN{t)]p 2 . (3.34) 

Although we have already formed joint stochastic equa- 
tions and averaged over unobserved processes, the above 
equations represent the step in Fig. |l| where we obtain 
the KSEs for classical probabilities of the APD states. 



6. Superoperator Kushner-Stratonovich Equation 



3.34 ) can now be 
-(3.16), using 



The expressions in Eqs. 
substituted into the RHS of Eqs. ( |3.14| ) 
Eq. ( |3.9|) for p. However, a simplifying deduction is that 
when the APD resets, the state of the quantum system 
must be equal to p 2 , which will already be normalized. 
We can therefore replace p by p 2 as the multiplier of the 
dM'it — Tdd) term. Similar reasoning allows us to ignore 
terms such as —dM{t)p 2 and j r Ti[pi]p 2 dt that will ap- 
pear in the equation for p 2 . The assumption of a known 
initial APD state allows us to make the following state- 
ment: at all times either Pq + Pi = 1 or P 2 = 1. After 
an avalanche, the detector is known to be in state 2 for 
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a time Tdd- Once the APD has been reset a realistic ob- 
server will not know if the APD is in state or 1 (but 
will know that it is not in state 2) until the time of the 
next avalanche. 

Also note that dN{t)dN'{t) = dN{t) = dN'(t) since 
here we are assuming that all avalanches are detected by 
the observer. This allows us to drop the prime on the 
non-Mar kovian resetting terms dN'(t — Tdd), rather than 
perform a repetitive conditioning on the observation of 
the resetting process, which is actually just inferred from 
the observation of an avalanche. 

Making these simplifications after substitution into 
Eqs. (p^-(fDl6|) gives 



dt (£ - 7dk - vJi c + M] + 7rTr[pi]) p 
- dN(t) Po + dN(t - r dd )p 2 (3.35) 
dt [(£ - 7 r + 7 r Tr[pi]) pi + {r/J[c + p] + 7dk) Po] 
-dM{t)p x (3.36) 

- dN '{t - r dd ) P2- (3.37) 



dp 
dpi 

dp 2 = dtCp 2 + dN (t) 



Pi 



Tr[pi] 



The above equations are the realistic quantum trajecto- 
ries based on the observation of avalanches in the APD. 
We term the equations SKSEs as we have obtained a 
quantum analog of the KSE in that from measurement 
we are conditioning the state of a supersystem that con- 
tains a quantum system, as explained in Sec. IC. The 
conditioned system state is 



PM = Po + Pi + Pi- 



(3.38) 



Note that Eqs. (|3~35|)-(^37|) are nonlinear, so that the 
state remains normalized. This is in contrast to the cor- 
responding equations in p4[ , which were linear and gen- 
erated unnormalized states. For simulation purposes, it 
is more convenient to use the unnormalized forms. 



C. Ignoring Imperfections 

Depending on the experimental situation, one or more 
of the imperfections in the APD may be negligible. Easy 
limits to take are those of rj = 1 and 7dk = 0, as these 
values can be substituted into the already derived equa- 
tions. More involved limits are those where 7" 1 and Tdd 
become small. These imply that one or two of the detec- 
tor states are superfluous. The first limit can be treated 
by a process of adiabatic elimination. 

When 7 r is much larger than the system rates (given 
by the eigenvalues of C) we can adiabatically eliminate 
state 1 of the dete ctor. When pi responds quickly to 
changes, p\ in Eq. ( |3.36| ) can be set to zero to give the 
slaved value 



Pi = —(j]J\c- 

7r 



p\ + 7dk)po- 



(3.39) 



Note that Cp\ and 7 r Tr[p 1 ]p 1 have been ignored as com- 
pared to 7 r pi and that Eq. ( 3.39| ) does not require the 



inclusion of a jump term due to avalanches (as these are 
inherited from the dependence upon p ). As p\ can now 
be determined from Eq. (3.39), the number of differen- 
tial equations for real variables that need to be evolved 
is greatly reduced 

If 7 r — * o o then state 1 can be ignored completely, with 
Eqs. (|3~35|)-(l3~37l) reducing to 



dp 
dpi 



dt (C - 7 dk - vJ[c + p}+ E[cW (t)]) po 
-dN(t)p + dN(t-T dd )p 2 , (3.40) 
dtCp2 — dN(t — T dd )p2 

(r)J[c + p} + 7dk)po 



- dN(t) 



Tr[(r? J[c + p] + 7dk) Po] 



(3.41) 



These equations describe the detector jumping straight 
from the ready state to the dead state. Thus, the exact 
time at which a charged e~-h + pair is created is known, 
but creation by a dark count or a photon is still indis- 
tinguishable. The statistics of dN (t) are now defined 
by E[dN(t)] = dtTr[(i]J[c + p] + 7dk)po] and we have 
PN = po + p 2 ■ Note that only one of po and p 2 will con- 
tribute to ps[ at any particular t ime. If the expression for 
F<[dN(t)] were included in Eq. (3.4C) then the —dt^ d ^po 
term would cancel due to it only having an effect on the 
transition rate out of state 0, which is now effectively 
being monitored. 

If, on the other hand, the dead time goes to zero, then 
state 2 of the detector can be removed. Thus, when 
an avalanche occurs, the detector jumps straight to the 
ready state. The relevant equations are then (with state 
1 not adiabatically eliminated) 



dp = dt (£ - 7 dk - r]J[c + //] + 7 r Tr[pi]) p 



dN{t)[p - 



Pi 



Tr[pi 



(3.42) 



dpi = dt[(C -7 r +7i-Tr[pi])pi + (r)J[c + p] + 7dk) Po] 
-dN(t) Pl . (3.43) 

Here, the stati stics of dN (t) are defined as for 
Eqs. ( 3.35 )— ( 3.37 ) but we now have p^ = po + pi- 

Finally, we consider 7 r — > oo together with Tdd = 0. 
This leaves an equation for p only. It is, 

dp = dt(C--y dk -r)J[c + p]+E[dN(t)])p + dN(t) 
(r)J[c + p}+ 7 dk) 



Ti[{r]J[c + p] +7dk)po^ 



- 1 



Po, 



(3.44) 



with E[dN 0Q1 defined as for Eqs. fl3.40| )-(|3^4T|). If 7dk = 
then Eq. ( 3.44 ) reduces to Eq. ( |2.17| ). Obviously, here 
we have p^ = Po- 



IV. PHOTORECEIVER 

When the incident photon flux is high, as in homodyne 
detection, an avalanche photodiodc cannot be used be- 
cause of the long dead time. Instead, a p-i-n photodiode 
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REALISTIC PHOTORECEIVER 



the presence of excess noise. We include only the Johnson 
noise Vj from the feedback resistor, which has a flat spec- 
trum 5,7 = Ak^TR. Here, fee is the Boltzmann constant 
and T is the temperature. This simplification (neglect- 
ing contributions from voltage noise of the operational 
amplifier and dark counts from the photodiode) can be 
justified for practical receivers with R ~ 10kf2 p4| . 

The output voltage V from the photoreceiver is given 
by sum of the filtered signal and the Johnson noise 



V = V + Vj. 



(4.3) 



FIG. 4: Homodyne detection by a realistic photoreceiver of a 
TLA placed in an optical cavity. The photocurrent generated 
by a below breakdown p-i-n junction, of efficiency 77, is input 
to a transimpedance amplifier. The voltage measured in the 
laboratory is V. 



is more appropriate. The inclusion of the high resistivity 
intrinsic region, i, allows the widening of the depletion 
region, which facilitates the absorption of light Q . Be- 
cause there is no avalanche, an external amplifier, such 
as a transimpedance amplifier (see Fig. must be used 
We consider the whole system of photodiode plus 
amplifier to constitute the photoreceiver. 

When a photon strikes the depletion region of the p-i-n 
junction, an electron-hole pair (that does not recombine) 
is produced with probability equal to the quantum effi- 
ciency 77. The charge carriers drift under the influence of 
the below-breakdown reverse bias, and the resultant cur- 
rent I is fed into an operational amplifier (op-amp) set 
up as a transimpedance amplifier. This has a low effec- 
tive input impedance, so that the diode acts as a current 
source, and / is converted into a voltage drop V across 
the feedback resistor, R. The capacitor C, in parallel 
with R, represents the total capacitance from the output 
of the op- amp back to its input, including capacitance 
added deliberately for the smoothing of noise and oscil- 
lations. If no electronic noise were present, the output 
voltage of the photoreceiver would be a filtered version 
of the input signal given, in the frequency domain, by 



-I{uj)R 
l + ituRC' 



(4.1) 



where the negative is included because the output of the 
op-amp will be such that the input is kept at virtual 
ground. 

It should be noted that if this were the case (that is, if 
there were no noise) then the input / could be perf ectly 
reconstructed from the filtered signal V, since Eq. (4.1) 
is equivalent to 



J(t) = -CV(t) - 



R 



(4.2) 



Thus the resultant quantum trajectories would be no dif- 
ferent from those of a photoreceiver with infinite band- 
width, since the input current to the amplifier is deter- 
minable. Everything of interest results therefore from 



Our aim is to find the quantum trajectory for the sys- 
tem, conditioned on continuously monitoring V. Since 
the voltage V, which describes the detector state, is a 
continuous variable, in this case S = M, the real line, and 
the supersystem can be described by an operator func- 
tion p(V). Finding the stochastic equation of motion for 
p(V) is quite involved. 



A. Realistic Quantum Trajectories 

In the following derivation we follow the order schemat- 
ically illustrated in Fig. [I] more closely than in the case of 
the APD. The averaging over the unobserved processes 
is performed as a final step. 

We begin by taking the output current I of the photo- 
diode to be that from a perfect (apart from its efficiency 
77) unbalanced homodyne detection of the output field of 
the system. For a LO tuned to the resonant frequency of 
the system u)q, of power P, and phase the current is 



I = e\JV /huj [ry(e 



— i<E> „ 



Vvm] , (4.4) 



where we have ignored the DC component due to the 
LO power and e is the electron charge. This current 
(in amps) is just the scaled current J(t) introduced in 
Eq. Q2.20D , multiplied by ey/P/Huj . 

The Gaussian white noise MM in Eq. (4.4) can be con- 



sidered to arise from two independent sources: the Pois- 
son statistics of the LO and the vacuum noise introduced 
by the inefficiency of the photodiode. The distinction 
between noise from these two sources will be illustra- 
tive when numeric simulations of the trajectories are per- 
formed (in a later paper). For now we need only consider 



1. Stochastic Master Equation 

The evolution of the quantum system condi tione d on 
I is given, in terms of the noise by Eq. ( [2.21 ) and 
is restated here for the reader's convenience 



d Pl = dt{C + y/v^W[e~ i9 c]}p 



(4.5) 
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The subscript J in Eq. (2.21) has been changed to /. The 
fact that p on the RHS of Eq. ( fL5| ) would usually be con- 
ditioned on previous values of I(t) is to be understood. 
A realistic observer does not have direct access to I{t), 
so later we will average over it. 



2. Stochastic Differential Chapman-Kolmogorov Equation 

N ow, Eq. ( |4.l| ) is equivalent to the stochastic equation 
(4.2). Since the voltage V is not directly measured, we 



must consider a distribution P(V) for it. Assuming tha t 
C > 0, and, for the moment, that I is known, Eq. ( |4.2| ) 
can be converted to an Ito stochastic Fokker-Planck equa- 
tion for P(V) conditioned on the photocurrent. This is 
done using the theory of Sect. II A, with the result being 



dPi(V) = 



d V + IR Prje 2 d 2 



dV RC 2huj C 2 dV 2 



P(V)dt. 



(4.6) 

A Fokker-Planck equation is one specific type of the 
broader class of di ffer ential Chapman-Kolmogorov equa- 
tions. As in Eq. (4.5) we are using the convention that 
subscripts indicate that the increment is conditioned on 
that result. For example, the probability distribution for 
V given a curr ent I is Pi(V) = P(V\I). It should be 
noted that Eq. (4.6) is different from a standard Fokker- 
Planck equation in that I still contains a white noise 
term. It is essential to reta in t his term as it is the same 
noise that appears in Eq. (4.4) - this will be used later 
in our derivation to link the evolution of the quantum 
system with that of the circuit of the photoreceiver. 



3. Kushner-Stratonovich Equation 

Next we need to determine the effect of the measure- 
ment of V on P(V). This can be calculated by using 
Bayes' theorem 



Pv(V) 



P(V) ' 



(4.7) 



Remember ing that the Johnson noise is white, it follows 
from Eq. (4.3) that Pv(V) is a Gaussian with mean V 
and variance 4k B TR/dt. This enables us to calculate the 
properties of P(V) from 

P(V) = J dVP v {V)P{V) (4.8) 
= (27T/3)- 1 / 2 / dVexp(-(V -V) 2 /2f3)P(V)(4.9) 



dVexp[(V- (V»(2V- V- (V))/20]P(V) 

x (2^/3)- 1 / 2 exp(-(V - (V)) 2 /2/3). (4.10) 

For convenience, the variance 4k B TR/dt in the above 
equations has been rendered by (3. For an expansion 



of the exponential i nside the integral to be possible the 
factorization of Eq. ( 4.10 ) had to be made. The variance 
in V is much smaller than that of V, which is of order 
0(1/ dt), due to the presence of the white Johnson noise. 

Expanding terms in the above expression to leading 
order, we find eventually that 



P(V) = 



exp(-(V-(V/)) 2 /2/3) 



1 + (eft 3 / 2 ) . (4.11 



It follows that we can write, correct to 0(1), 

■dWj(i) 



V = (V)+ y/Ak B TR- 



dt 



(4.12) 



where dWj(t)/dt is another Gaussian white noise source, 
independent of Despite the fact that it arises from 
the Johnson noise, it is not the actual Johnson noise. 
It represents noise that the observer can determine from 
the measurement result and (V). Labeling the actual 
Johnson noise by dWj(t) we have 



V = V + y/4k B TR 



dWj(t) 
dt ' 



so that 



dWj(t) = dWj(t) + dt 



(v-(v)) 

V4fc B Ti? ' 



(4.13) 



(4.14) 



Substitution of Py(V) and P(V) into Eq. (|4.7|) gives to 
first order in dt 



P V (V) = P(V) + dWj(t) V { l ] P(V). 

V4kbJ R 



(4.15) 



That is, we have derived the Kushner-Stratonovich equa- 
tion 

dP v (V) = dW,(t) {V _J^ P(V), (4.16) 

where the conditioning upon V is contained in dWj (t) = 
dt (V - (V)) /y/ik B TR. This quantity is often referred 
to as the residual or innovation [[l8], |l9| . 

4- Joint Stochastic Equations 

To see how V conditions the quantum system, we form 
the quantity p(V) = pP(V), where p is here independent 
of P(V) because we are still imagining / to be known at 
all times. The time evolution of p(V), given that V and 
/ are known, is found from 

p{V)+dp iy {V) = [P(V)+dP!(V) + dP v (V)] 

x(p + d Pl ), (4.17) 



with the use of Eq. (4.6), Eq. (4T9) and Eq. (O) 
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Substitution into Eq. (4.17) gives 
iV-{V)) 



Piy{V) 



l + dWj(t) 



dt 



d V + IR Prje 2 



d 2 



<>}- nc - < 'V 2 ') 



P(V) 



x [l + dt(C + ^(t)n[e-^c])]p. (4.18) 

The reader is reminded that is uncorrelated with 
dWj(i). Upon using / given in Eq. (i.4) and noting 
that it co ntains the white noise £(i), the expansion of 
Eq. dy||) to 0(dt) is 



+ {^ew + [^(*)] 9 ?V^} 

x [e~^cp{V) + e^p{V)c^ - ( e -**c + e*V) p(V)] , (4.19) 



where the form of the superoperator Ti in Eq. ( 2.19 ) has 
been used. 



5. Average Over Unobserved Processes 

Finally, in reality V is known but / is not. Therefore 
we should average over the noise £(t), but keep the volt- 
age noise dWj(t). The relevant averages are Eq. (2.2), 
implying di 2 E[£(t)£(t)] = dt. An important consequence 
(necessary for consistency) is a cancellation of the terms 



involving (e 



). For convenience, we define a 



dimensionless voltage v = V \J C/AIcbT, a rate 7 = 1/RC 
and a dimensionless noise power 



N = 



AkBTtkuo 
r]RPe 2 



(4.20) 



This last expression is the ratio of the Johnson noise 
power {^UbTR) to the low- frequency power in V due 
the noise in the photocurrent {R 2 e 2 V-q/hw Q ). The lat- 
ter noise is that which would be present if a vacuum field 
was being combined with the LO field at the homodyne 
detection beam splitter. 

Averaging over £(i) and using v, 7 and N we then 
obtain the following superoperator Kushner-Stratonovich 
equation (SKSE) for p(v): 



dp v (v) 



dv V N 



<-y (ft Q 

2Nd^ 2+1 d-v V]p[v) 



cp(v) + e l ® p{v)J] 




v))p{v). 



(4.21) 



As noted below Eq. ( 4. 16| ) , the conditioning on the quan- 
tity that is measured in the laboratory, V, is contained 



in dWj(t), with 



^dW J (t)=dt 1 (^^V-(v) 



(4.22) 



Although the description of the state of the photore- 
ceiver (by v) is essential in obtaining realistic trajecto- 
ries, it is the evolution of the quantum system that we 
are most interested in. This state is given by 



Pv 



p(v)dv. 



(4.23) 



The average (v), which appears in Eq. ( 4.21 ), is found 
from 



v) = dv Tr[p(v)]v. 



(4.24) 



B. Effective Bandwidth 



As mentioned earlier, if there were no excess (Johnson) 
noise added to the input signal, then this input signal 
could be perfectly reconstructed from the filtered out- 
put. However, in the case of realistic detection presented 
above there is noise, characterized by N, linearly added 
to the output of the filter. The extent to which informa- 
tion is lost due to this noise depends on the magnitude 
of the noise, the filter bandwidth, and the nature of the 
evolution of the monitored system. In this subsection, 
we argue that in the limit of small noise N <C 1, the 
quality of the photoreceiver can be characterized by an 
effective bandwidth that depends upon the noise and the 
filter bandwidth. We expect the information loss to be 
small only if this effective bandwidth is large compared 
to the relevant system frequencies. 
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In terms of physical parameters 

e / 

B= 2Cyk B TRh^ - (4 - 2?) 

The photon flux P / Hujo of the LO cannot be increased ar- 
bitrarily as the regime in which the photodiode responds 
linearly to the field incident upon it must be adhered to. 
The resistor R cannot be reduced to zero as it is neces- 
sary for the reduction of the effective input impedance of 
the op-amp so that the photodiode can act as a current 
source. It also converts the current to a voltage, which 
presumably would be swamped by other, not considered, 
noise sources if reduced too far. Minimizing stray capaci- 
tances and the temperature of operation will improve the 
photoreceiver operation. 



FIG. 5: The dimensionless quantities J 2 /(j 2 + lo 2 ) and TV are 
plotted, showing their frequency dependence (or lack of it). 
In this plot 7 = 1.5 and TV = 0.1. The frequency at which the 
vacuum signal drops in power to that of the noise is equal to 
the effective bandwidth B. 



We identify the effective bandwidth, B, as being 
roughly the frequency at which the noise in a vacuum 
input signal is 'lost' in the Johnson noise. That is, the 
frequency at which the power of the vacuum input signal 
becomes equal to the Johnson noise power (which has a 
flat spectrum) . Obviously the loss of very high frequency 
noise will invariably occur in any practical photoreceiver. 
The question of whether this loss is important or not is 
answered by looking at the eigenvalues of the Liouvillian. 
If the input noise is lost at frequencies that are well above 
the system rates then this noise would not significantly 
affect the evolution of the system. However, if noise is 
lost at frequencies at which the system can respond to, 
then the purity of the system state will decrease consid- 
erably. This is why we expect the effective bandwidth, 
as defined above, to be the relevant parameter for the 
photoreceiver. 

Equating the vacuum signal and Johnson noise powers 
at the frequency u> — B gives 



7 



1 2 + B 2 



Ak B TRHuj Q 
R 2 e 2 Prj 



TV 



(4.25) 



This is shown graphically in Fig. |5|. Obviously if TV is 
too large then this equation has no solution in the real 
numbers, and indeed the intuitive picture behind the ar- 
gument fails. But for small TV we have 



B = 



7^1 - TV 7 



TV 



TV 



(4.26) 



This result will be investigated further for two different 
quantum systems in the following paper. 



C. Adding Noise Only 

The limit in which the power of the Johnson noise goes 
to zero has already been discussed, with the result being 
perfect detection. If the response time of the RC circuit 
vanishes then, provided Johnson noise is still added, the 
input signal is still obscured. For this to be a physically 
sensible consideration we assume that the capacitance, 
rather than the resistance, is zero. 

The output voltage of the photoreceiver is then given 

by 



V 



-IR + \/4k B TR 



dWj(t) 
dt ' 



(4.28) 



as indicated by the form of Eq. (4.3). Here dWj(t)/dt 
represents the 'real' Johnson nois e whi ch is of a Gaussian 
nature. Using Eq. (Q|) and Eq. ( 4.20 ) gives (with $ = 
for convenience), 



V 



n 



1 + TV 



(c + c f ) + 



V 



1 + TV 



0(t), (4.29) 



where 9(t) is normalized white noise consisting of a linear 
combination of the vacuum in put noise and the Johnson 
noise. Comparison with Eq. (4.4) shows that the effect 
of adding noise but no filtering is to reduce the efficiency 
by a factor of 1 + TV. It can be shown using Bayesian 
analysis that this result holds for the form of the quantum 
trajectory as well. 



V. CONCLUSION 

The evolution of open quantum systems conditional 
upon detection results from realistic detectors cannot in 
general be generated by standard quantum trajectory 
theory (stochastic master equations). A method for gen- 
erating this evolution was proposed in Ref. 1 14 . The cru- 
cial element is to include correlations between the system 
and classical detector states which cannot be observed in 
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practice. In this paper we have given a full description of 
this method, and shown how it can be applied in quan- 
tum optics. In particular, we derive realistic quantum 
trajectories for conditioning upon photon counting using 
an avalanche photodiode, and homodyne detection us- 
ing a photoreceiver. These equations were presented in 
Ref. |l4j, but a full derivation was not given there. 

In this paper we have not provided any solutions of 
the equations we have derived. To find and study these 
solutions is not a minor task, which is why it is reserved 
for the following paper p3| ]. There we use our equations 
to determine the evolution of a two-level cavity QED sys- 
tem, conditioned on four different types of detection (us- 
ing the two detectors mentioned above). We also consider 
another system which can be treated analytically for re- 
alistic homodyne detection. Our study achieves five im- 
portant aims, which are outlined in the introduction to 
that paper. 

It is worth re-emphasizing the generality of our ap- 
proach. It is applicable not just in quantum optics, 



but in all areas where open quantum system theory is 
used. Prime examples are the mesoscopic measurement 
devices used in quantum electronics, such as single elec- 
tron transistors. This device has received much atten- 
tion as a possible measurement device for solid state 
qubits |55|, |3(| [57). Conditional states may be used in 
quantum computation both for preparation, and for non- 
deterministic gate implementation |38). In this context, 
the necessity of being able to relate realistically available 
results to the state of the qubit is obvious. 

The greatest significance of our work is in the field 
of quantum control. The conditional state of a quan- 
tum system is synonymous with an observer's knowledge 
about that system. Thus (assuming all uncertainties are 
properly taken into account) it is by definition the opti- 
mum mathematical object to use to control that system. 
Properly taking into account detector imperfections is es- 
sential to building optimal control loops. We thus expect 
our theory to have broad applications in future quantum 
technology. 
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